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ABSTRACT 

Recent analyses of multiwavelength light curves of gamma-ray bursts afterglows point to 
values of the magnetic turbulence well below the canonical ~ 1 % of equipartition, in agree¬ 
ment with theoretical expectations of a micro-turbulence generated in the shock precursor, 
which then decays downstream of the shock front through collisionless damping. As a direct 
consequence, the Compton parameter Y can take large values in the blast. In the presence 
of decaying micro-turbulence and/or as a result of the Klein-Nishina suppression of inverse 
Compton cooling, the Y parameter carries a non-trivial dependence on the electron Lorentz 
factor, which modifies the spectral shape of the synchrotron and inverse Compton compo¬ 
nents. This paper provides detailed calculations of this synchrotron-self-Compton spectrum in 
this large Y regime, accounting for the possibility of decaying micro-turbulence. It calculates 
the expected temporal and spectral indices a and customarily dehned by Fy oc 
various spectral domains. This paper also makes predictions for the very high energy photon 
flux; in particular, it shows that the large Y regime would imply a detection rate of gamma-ray 
bursts at > lOGeV several times larger than currently anticipated. 
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1 INTRODUCTION 

The afterglow radiation of gamma-ray bursts, which spans the 
range from the radio domain to the X-ray and possible higher fre¬ 
quencies, with a characteristic decrease in time of the peak flux, 
is nicely interpreted in terms of the synchrotron radiation of ultra- 
relativistic electron that have been accelerated at the forward ultra- 
relativistic collisionless shock wave of the outflow, see e.g. Piran 
(2004) for a detailed review. 

This model thus opens a beautiful connection between as¬ 
tronomical data and the microphysics of the rather extreme phe¬ 
nomena that relativistic collisionless shocks represent. Zooming 
in on microphysical scales reveals the shock front as a micro- 
turbulent magnetic barrier which isotropizes the incoming back¬ 
ground plasma, see e.g. Moiseev & Sagdeev (1963) for pioneering 
studies and Kato & Takabe (2008), Spitkovsky (2008a) for detailed 
numerical simulations. This micro-turbulent barrier itself results 
from the build-up of electromagnetic micro-instabilities ahead of 
the shock front, in the shock precursor where supra-thermal parti¬ 
cles mix with the incoming background plasma, e.g. Medvedev & 
Loeb (1999), Wiersma & Achterberg (2004), Lyubarksy & Eichler 
(2006), Achterberg & Wiersma (2007), Achterberg et al. (2007), 
Bret et al. (2010), Lemoine & Pelletier (2010, 2011), Rabinak et al. 
(2011), Shaisultanov et al. (2012), Lemoine et al. (2014a, b). 

The resulting micro-turbulence, with typical length scale 
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Aas ~ c/tJpi (ojpi the ion plasma frequency) and typical strengtl0 
es ~ 0.01, thus emerges as a key ingredient for the micro¬ 
physics of collisionless shocks. Actually, it also represents a key 
requisite for the relativistic Fermi process, since this latter can take 
place (in ideal conditions) only when an intense micro-turbulence, 
with power on scales smaller than the gyroradius of the acceler¬ 
ated particles, is able to unlock the particles off the background 
magnetic field lines, by scattering them faster than a gyration 
time in this background field (Lemoine et al. 2006, Niemiec et 
al. 2006, Pelletier et al. 2009). In terms of magnetization g = 
j (Stt dT^nmpC^), with B the background field expressed in 
the downstream frame, this condition amounts to ct ^ e% 
(Lemoine & Pelletier 2010, 2011, Lemoine et al. 2014a), a con¬ 
dition which is indeed satisfied for the forward shock of gamma- 
ray bursts, since a ~ 10“® in the interstellar medium. This point 
of view has been confirmed by particle-in-cell simulations (e.g. 
Spitkovsky 2008b, Martins et al. 2009, Nishikawa et al. 2009, 
Sironi & Spitkovsky 2009, 2011, Sironi et al. 2013, Haugbplle 
2011 ). 

Finally, these micro-instabilities may also build the magne¬ 
tized turbulence in which the electrons eventually produce the af¬ 
terglow radiation in a synchrotron-like process (Medvedev & Loeb 
1999) - jitter effects are expected to be weak in the conditions 


^ eg = (577 ^/(Stt 4r^nmpC^) represents the magnetic energy fraction 
of equipartition, if Ft, 1 represents the relative Lorentz factor between 
upstream and downstream, n the upstream proper density. 
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typical of those shock waves - provided it survives collisonless 
damping downstream of the shock front (Gruzinov & Waxman 
1999). Recent analyses of the damping of this Weibel-type micro¬ 
turbulence yield a damping rate oc fc®, where k denotes a turbulent 
wavenumber, indicating that small scales are dissipated first, but 
that large scales may survive longer (Chang et al. 2008, Lemoine 
2015). In turn, this implies that the turbulence strength, or es, 
should decay as a power-law in (proper) time (or distance) down¬ 
stream of the shock front, with an index which itself depends on 
the (unknown) micro-turbulent power spectrum at the shock front. 
Interestingly, this time dependence of ts turns out to be encoded in 
the multiwavelength synchrotron spectrum, since electrons of dif¬ 
ferent Lorentz factors cool at different times, hence in regions of 
different magnetic field strengths (Rossi & Rees 2003, Derishev 
2007, Lemoine 2013). 

Although one cannot exclude that other external instabilities 
would pollute the blast with magnetized turbulence, it is tempt¬ 
ing to consider that the generation of micro-turbulence could be 
responsible at the same time for the formation of the shock, for 
the acceleration of particles and for the radiation of these parti¬ 
cles. It is furthermore tempting to follow this thread to use the 
multiwavelength spectrum of gamma-ray bursts as a tomograph of 
the magnetized turbulence. As a matter of fact, the recent detec¬ 
tions of extended emission of gamma-ray bursts in the > 100 MeV 
band by the Fermi-LAT instrument do point to a net decay of the 
micro-turbulence behind the shock front (Lemoine et al. 2013), 
with e_B oc (f6L)pi/100)“* and —0.5 < crt < —0.4. This argu¬ 
ment can be recapped as follows: if the accelerated particles scatter 
in a micro-turbulence, the maximal synchrotron photon energy is 
limited to a few GeV at an observer time of 100 s (Kirk & Re- 
ville 2010, Plotnikov et al. 2013, Wang et al. 2013, Sironi et al. 
2013); this maximal energy scales as the square root of eB+, i.e. 
the magnitude of the turbulence in the vicinity of the shock front, 
and the above value assumes eB+ ~ 0.01; therefore, the inter¬ 
pretation of this extended > 100 MeV emission as a synchrotron 
process points to the existence of a strong micro-turbulence close 
to the shock front; on the other hand, multiwavelength fits of the af¬ 
terglows for these Fermi-LAT gamma-ray bursts indicate that low- 
energy photons are produced in regions of rather low es_, of the 
order of 10“® — 10“®; this discrepancy between es_ and e_B+ is 
naturally interpreted as the decay of the micro-turbulence through 
collisionless damping. As discussed in Lemoine et al. (2013), the 
decay rate at ~ —0.5 to — 0.4 further matches the results of a 
detailed numerical experiment reported in Keshet et al. (2009). 

The decay of micro-turbulence has also been proposed as 
a possible solution to the abnormal spectral indices observed in 
the prompt emission phase (Derishev 2007), but admittedly, the 
physics of mildly relativistic shock waves may well differ from that 
of ultra-relativistic shock waves; in particular, the extended precur¬ 
sor size in mildly relativistic shocks opens the way to other insta¬ 
bilities operating on larger length scales. Although the present con¬ 
siderations can also be generalized to the case of internal shocks, 
all of the discussion that follows will focus on the external ultra- 
relativistic shock front. 

An interesting consequence of a decaying micro-turbulence, 
or more generally, of a low average value of eb as measured in the 
Fermi-LAT and other bursts (e.g. Kumar & Bamiol-Duran 2009, 
2010, He et al. 2011, Santana et al. 2014, Barniol-Duran 2014), is 
a large Compton parameter Y, at least if one omits the influence 
of Klein-Nishina effects, which actually depend on electron en¬ 
ergy hence on observed frequency, see below. Recall indeed that 
in the standard model, assuming that electrons cool through in¬ 


verse Compton interactions on their synchrotron spectrum in the 
Thomson limit (i.e. neglecting Klein-Nishina - KN - effects), 
Y ~ {rjeeltBY^^ SilY > 1, with ~ min [l, ( 7 c/ 7 m)^“^] 
the cooling efficiency of electrons, 7 c denoting the Lorentz factor 
of electrons which cool on a dynamical timescale, 7 in the minimum 
Lorentz factor of the injected electron power-law and —p the index 
of this power-law, e.g. Panaitescu & Kumar (2000), Sari & Esin 
(2001), Piran (2004). This Y parameter also reflects the power of 
the inverse Compton component relatively to the synchrotron com¬ 
ponent, therefore the ratio of emission at very high energies to that 
in the X-ray range. The possibility of a large Y parameter should 
thus increase the chances of observating gamma-ray bursts at very 
high energies with upcoming Cerenkov telescopes. The observa¬ 
tion of this inverse Compton component would then open a new 
spectral domain with which one could study the microphysics of 
the turbulence in the blast. 


However, the computation of the synchrotron-self Compton 
(SSC) spectrum of the blast in the large Y regime is not triv¬ 
ial because the shape of the synchrotron spectrum influences the 
cooling history of the electrons, which itself determines the syn¬ 
chrotron spectral shape. In particular, the KN suppression of the 
inverse Compton cross-section may itself modify this cooling his¬ 
tory, hence modify the very shape of the synchrotron spectrum, see 
Nakar et al. (2009), Wang et al. (2010); see also Barniol Duran & 
Kumar (2011) for the particular case of GRB090902B and Bosjnak 
et al. (2009), Daigne et al. (2011) for related discussions in the con¬ 
text of gamma-ray burst prompt emission. In the case of decaying 
micro-turbulence, this issue is more acute, because the synchrotron 
power vFi, may be rising with frequency, due to the fact that lower 
frequency photons are emitted by electrons of longer cooling time, 
in regions of lower magnetic field strength (Lemoine 2013). There¬ 
fore, the Lorentz factor-dependent KN limit determines the radia¬ 
tion intensity on which an electron can cool. Finally, so far the in¬ 
fluence of decaying micro-turbulence on the synchrotron spectrum 
has been studied for a fixed Y parameter, independent of electron 
Lorentz factor, hence the general shape of the SSC spectrum is not 
known in this physically relevant case. 


These considerations motivate the present study, which calcu¬ 
lates the SSC spectrum for a relativistic blast wave in the large Y 
regime. Although the prime objective is to understand how the de¬ 
caying micro-turbulence affects the SSC spectrum, the present dis¬ 
cussion also addresses the case of a uniform, low value of eb. The 
results are applied to the case of the GRB afterglows, by calculating 
the spectrum at various observer time and by plotting, in particular, 
the spectral slopes in various frequency windows. The spectrum is 
also evaluated in the multi-GeV range, with a proper account of 
synchrotron and SSC contributions, to make clear predictions for 
future Cerenkov observatories such as HAWK and CTA. 


This paper is organized as follows: Sec. [^presents an analyti¬ 
cal description of the spectrum in the slow-cooling limit, account¬ 
ing for KN effects, and introduces a fast algorithm to compute this 
spectrum in the fast cooling regime; Sec.j^then computes the light 
curves and plots the temporal and spectral indices a and /3, com¬ 
monly defined by Fj, oc and gives predictions for the 

very high energy photon flux; finally. Sec. [^summarizes these re¬ 
sults and provides some conclusions. 
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2 SSC SPECTRUM 
2.1 Set-up 


The set-up is as follows: electrons are swept-up by a relativistic 
shock front propagating in a density profile n oc then acceler¬ 
ated on a short time-scale to a power-law dA^e,o/d7 c< 7“^0(7 — 

7m) above a minimal Lorentz factor ym- The minimal Lorentz fac¬ 
tor is defined as usual by 7m = tePhmp/me{p — 2)/(p — 1), 
with Ce ~ 0.1. Cooling takes place on much longer timescales. 
When describing the cooling history, the total inverse Compton 
cross-section is modelled as a top-hat, with a = ut for u < 
and zero otherwise, vipf) denoting the frequency of photons with 
which electrons of Lorentz factor 7 interact at the KN limit (Nakar 
et al. 2009, Wang et al. 2010), i.e. 


FbmeC^ 

/i 7 (l -I- z) 


( 1 ) 


All frequencies are written in the observer rest frame; as mentioned 
above. Lb denotes the Lorentz factor of the blast in the source rest 
frame. The Compton parameter can then be approximated as 


where 


y( 7 ) 


t/rad [V < ^'(7)] 


t/rad [< £'(7)] = / 

Jo 


( 2 ) 

(3) 


represents the comoving radiation energy density at frequencies 
V < £'(7), on which the electron of Lorentz factor 7 can cool. Of 
course, if iz > tZpeak at which the differential energy density i7i/(zz) 
reaches its maximum, then Urad [v < £'( 7 )] ~ ~ 

t/rad and it does not depend on 7 (t/rad denotes here the total co¬ 
moving energy density). In principle, t/rad includes all form of ra¬ 
diation, synchrotron and inverse Compton alike; however, multi¬ 
ple Compton scattering can be neglected for standard gamma-ray 
bursts parameters, hence t/rad is hereafter determined with the syn¬ 
chrotron spectrum only. 

The quantity Ub{Y represents the energy density contained 
in the magnetic field and it depends on 7 if the turbulence de¬ 
cays in (proper) time behind the shock, because electrons of dif¬ 
ferent Lorentz factors then cool in regions of different magnetic 
field strengths. In this work, the decay law of the turbulence takes 
the power-law form esit) oc far from the shock front and 
es ~ es+ = 0.01 close to the shock front. How far is ex¬ 
pressed in terms of the proper time t since the plasma element was 
injected through the shock, i.e. t = x//3d in terms of the down¬ 
stream (comoving) distance to the shock front x and /3d, the shock 
front velocity relative to the downstream rest frame. According to 
PIC simulations, the characteristic (temporal) scale A separating 
far from close to the shock front is A ~ 10^“® cUpi, see Chang et 
al. (2008), Keshet et al. (2009) for simulations, as well as Lemoine 
(2013), Lemoine (2015) for discussions of this issue; for reference, 
~ 7.5 x s for a relativistic blast wave propa¬ 

gating in a medium of proton density nocm~®. The uncertainty 
on the value of A does not really influence the results presented 
below, because it can be embedded in that associated to the de¬ 
cay power-law exponent at - Recent work mentioned above suggest 
Q.t ~ —0.5 —>■ —0.4 for a decay law es ~ e_B+ (ta;pi/100)“*, 
hence the following adopts A = 100cUp3 and at = —0.4. This 
decay implies that the minimum magnetic field in the blast, close to 
the contact discontinuity, is characterized by the equipartition frac¬ 
tion es_ = es+ (tdyn/A)“* in terms of the dynamical timescale 


fdyn = r/(Lbc). Typical values are tdyn ~ 10® — 10® s at an ob¬ 
server time of 10"^ s (e.g. Fb ~ 20 —30 and r ~ 10^^ —10^® cm), 
with a mild dependence on the model parameters, leading to val¬ 
ues of the order of es- ~ 10~® for at ~ —0.5. Note that 
tdyn/A oc evolves slowly as a function of observer 

time, hence so does es_. 

In this work, it is assumed that a particle emits its synchrotron 
radiation at the location at which it cools, if it cools on a dynam¬ 
ical timescale; in the opposite limit, it is assumed that it emits its 
synchrotron photons in the magnetic field close to the contact dis¬ 
continuity, of strength SB- (associated to the parameter es_). The 
justification of this approximation is as follows: the energy emitted 
in synchrotron photons by a particle up to time t, along its trajec¬ 
tory downstream of the shock, can be written as 

J dr5B^(r)7e(T)^e(r) (4) 

as a function of the time dependent Lorentz factor 'yei'r) of the 
particle along its cooling trajectory, with fieir) ~ 1 the particle 
velocity in units of c. One can then show that Esyn is dominated 
by the contribution at i ~ £cooi(7), where 7 denotes the initial 
value 7e(0), as follows. At early times t <C icooi(7), 'yei'’') ~ 
7 and the integrand scales as Ub{t) oc r“*, hence the integral 
is dominated by the large time behavior if at > —1, which is 
indeed an explicit assumption of the present work. At late times 
T ^ /cool ( 7 ), one can obtain an upper bound on how fast 7e('r) 
decreases by considering synchrotron losses only, i.e. by neglecting 
inverse Compton losses. The integration of 


d7e 

dr 


4 

3 


crTc[/s(r) 

rrisC^ 


7e(^)^e(r)^ 


(5) 


leads to 7e(r) ~ 7 [r//cooiM] ^ for r > /cool (7)- 

Therefore, the integrand in Eq. (H behaves as for r > 

/cooi( 7 ), hence Eayn is indeed dominated by the contribution at 
/ ~ /cool (7) provided Qt > -1. 

This thus supports the above approximation that a particle 
with initial Lorentz factor 7 > 7c emits its synchrotron radia¬ 
tion on a magnetic field of strength SB [/cool (7)], with /cool (7) the 
cooling time of the particle. Of course, if /cool (7) > /dyn, the 
particle does not actually cool, and Esyn is then dominated by the 
upper bound / ~ /dyn, i.e. particles radiate synchrotron radiation 
in the relaxed magnetic field at the back of the blast, of strength 
SB-. 

Since /cooi(7c) = /dyn by definition, and since ( 5 / 3 (/dyn) = 
SB -, also by definition. 


Us (7) ~ Us-max 1, 


^ /cool (7) 
V /dyn 


Of 


( 6 ) 


with Us- = SB^ /{Stv) = Us ( 7 c). 

In line with the above discussion, electrons of Lorentz factor 
7 radiate their energy through synchrotron at a typical frequency 


t^syn( 7 ) « [/cool (7)1 7 ^ - t^c 

(7) 

Here, Sc = Ssyn (7c) denotes the synchrotron peak frequency for 
particles of Lorentz factor 7 c. 

Finally, the cooling timescale can be written 


/cool (7) 

/dyn 


t/2 


/cool (7) 


/d 


/d 


1 + Ys Ub- 7 c 
i + y(7)UBWT ’ 

i-f y(7)1"^/^^+“*’ 


l-FUc 



-l/(l + Ot) 

(8) 
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The second equality is obtained by replacing ?7s ( 7 ) jUs- with its 
value given in Eq. assuming 7 > 7 c. 

Equations 0 and 0 then allow to derive the follow¬ 

ing scalings for 7 > 7 c: 


(7b ( 7 ) 

■i + y(7) 

Ub- 

L i + yc 

'=yn(7) 

i + y( 7 )l 

Uc 

i-byc 




— OLt/[2{l + OLt)] 




-oct/(l-\-oct) 


(9) 


2 -at/[2{l + c«t)] 

( 10 ) 


yc(l + yc)-“‘/(l + “‘) “ C/rad(t2<Z2c) Uc/ 

( 11 ) 

Yc stands for Y (jc)- 

Once the ratio ( 7 rad(t 2 < v)/U- cz.a(v < Vc) - which captures 
KN effects - has been specified, it is possible to derive the scalings 
of Y ( 7 ) as a function of 7 , hence of r'syn( 7 ). One should empha¬ 
size that in the above expressions, 7 represents the initial Lorentz 
factor of the particle, after acceleration has shaped the power-law, 
but before cooling has effectively taken place. 


2.2 Synchrotron spectrum 


Standard calculations of the synchrotron spectrum of a blast wave 
generally derive the stationary electron distribution in the blast, by 
solving a transport equation in momentum space, averaged over the 
depth of the blast, accounting for injection of the power-law at the 
shock and for cooling downstream (e.g. Sari et al. 1998, Sari & Esin 
2001 ); this yields a standard broken power-law shape dNs/dy^ 
with indices —2 for 7 c < 7 < 7 m or —p for 7 m < 7 < 7 c 
and —p — 1 for max( 7 c, 7 m) < 7 . However, one can also de¬ 
rive the synchrotron spectrum by a direct mapping of the energy 
initially stored into the electron population to that radiated in syn¬ 
chrotron; this approach is simpler in the present case and it works 
as follows. 

Electrons injected with Lorentz factor 7 > max( 7 m, 7 c) 
emit a fraction [1 -|- ^( 7 )]”^ of their energy in synchrotron ra¬ 
diation; the energy density contained in such electrons (within 
an interval din 7 ) is itself a fraction (p — 2) ( 7 / 7 m)^~^ din 7 
of the energy density Ue contained in electrons immediately be¬ 
hind the shock front, before cooling has started to take place; 
Ue = Ce YT^nnipC^. Consequently, the synchrotron flux received 
at frequency v = r'syii( 7 ) can be written: 

jp 1 4^2 47 rr^c(p - 2 )f/e f 7 7 ”^ dln 7 

- 4^3 l + y(7) V7/ dW 

( 12 ) 

The above assumes p > 2 , but it can be generalized to any in¬ 
dex p > 1 ; indeed, the above picture assumes that particles with 
7 > max( 7 m, 7 c) radiate their energy at I'syni'y) then do not con¬ 
tribute anymore to the synchrotron spectrum, which is a reasonable 
approximation if p > 1 ; in contrast, if p < 1 , the radiation of the 
high energy particles during their complete cooling history domi¬ 
nates that of the lower energy ones. 

Similarly, one can write down the flux for 7 m > 7 > 7 c or 
7 c > 7 > 7 m, whichever occurs, as follows: 


syn — r^mf^i^nnsyn 


7 1 + 


1 + Y 

(3-p)/2 


(7m > 7 > 7c) 

(7c > 7 > 7m) 
(13) 


with Ym = Y(jm) and i/m = i^synijm)- The spectrum at 
< tZc is indeed unchanged with respect to the standard case, 
although the value of Ub- must be used to compute the character¬ 
istic frequencies. Regarding the fast cooling limit 7 c < 7 < 7 m, 
the factor can be understood by noting that all injected particles 
with 7 > 7 m shift from 7 to 7 c during their cooling history, and 
that at each point along this cooling trajectory they radiate (in syn¬ 
chrotron) a fraction ( 7 / 7 m) (1 -b Ym) / (1 -b y) of the energy ra¬ 
diated (in synchrotron) by a particle of Lorentz factor 7 m. Note 
also that in the limits at —>■ 0 , y( 7 ) <C 1 , one recovers the usual 
scaling i'Fi,^syii oc since v oc 7 ^. 

Provided 7 c and yc are given, one can derive r'ayn( 7 ) and 
y( 7 ), hence 'y{i') and Y{u) by using Eqs. 0 , and ([^, 
then vF^^syn using Eqs. l |12| l and ( |13| l. Consider as an example 
the case Y ^ 1, omitting KN effects, i.e. v —> -boo for all 
u: for 7 > 7 c, one derives y( 7 ) oc 7 “* from Eq. jl 1[ (, hence 
Up oc 7 ^““*'^^ from Eq. 1 10 1 , hence i/Tl,,syn oc 7 ^“^/(l + y) oc 
^{ 2 -p-at)i( 2 -at/ 2 ) ^ j 21 . xhis latter scaling matches the 

explicit calculation of Lemoine (2013), which computes the inte¬ 
grated cooling history of the electron power-law in the decaying 
magnetic field. This example also confirms that uF,^^syii rises with 
u if at < 2 — p, hence KN effects should not be omitted in an 
accurate calculation of i 2 T 7 ,syn. 

The above provides the tools needed to calculate the syn¬ 
chrotron spectrum; explicit calculations in both the slow and the 
fast cooling regimes are provided in the following sections. 


2.3 Inverse Compton spectrum 

In order to calculate the inverse Compton emissivity, the syn¬ 
chrotron spectral density Tj/,syn must be folded over the Compton 
cross-section and the cooled particle distribution, integrated over 
the blast. The following derivation of the cooled distribution func¬ 
tion relies on the observations that dfcooi ( 7 ) /dx < 0 and that the 
cooling history 7 cooi(f) at times t > to does not depend on the 
history at times t < to, i.e. electrons of various Lorentz factor 
follow a universal cooling trajectory 7 cooi(f). 

One writes dA^e.o the (comoving) rate at which electrons are 
swept-up and accelerated into a power-law in a Lorentz factor in¬ 
terval dy. One also defines dA^e, which represents the total number 
of electrons in the blast in a Lorentz factor interval d 7 e; note the 
difference of notation: 7 e refers to a value of the Lorentz factor of 
the blast averaged distribution, while 7 corresponds to the Lorentz 
factor of the injection distribution, i.e. before cooling has occurred. 

The injection of particles with Lorentz factor 7 populates the 
downstream with particles of Lorentz factor 7 e = 7 over a frac¬ 
tion fcooi/fdyn of the depth of the blast. In contrast, particles in¬ 
jected with 7 c > 7 > 7 m retain their Lorentz factor (no cooling) 
and populate the whole blast. Therefore, one can write the average 
distribution for 7 e > 7 m as: 


dA^e 

d7e 


dA^e, 

d7 


'(7 


7e) min [fdyn, fcool(7)] (7 > 7m) . 


(14) 

Note also that fdyn sets the scale for adiabatic losses, while fcooi 
sets the scale for radiative losses. 

In the fast cooling regime, all electrons injected with Lorentz 
factor 7 > 7 m cool down to 7 c, following a same cooling history 
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7 cooi(t), therefore 


dN, 

d7e 




dA^e, 


d7 


be 


fcool( 7 e) 
^^e,0 - 

7e 


7cooi('r)] 


(15) 


where o = djdNe o/d 7 represents the total injection 

’ bm ’ 

rate of electrons. The last expression in Eq. jl5[ l follows from the 
definition fcool( 7 e) = 7 e |d 7 e/df|~^ after interchanging the inte¬ 
grals. 

Using the scaling of fcooi( 7 ), one then derives 


diVe 

d7e 

diVe 

d7e 


dNe 

dje 


fVe.ofdyn / Ji 


7n 


7n 


A^e,ofdyn ^ 7e 
V7n: 


7ix 


-1 


1 + r(je) 


Jl 

7. 


(7m < 7e < 7c) 

— 1/ (l+at) 


l+Uc 

fVe,ofdyn f 7^ 


7rx 


7u 


1 + y(7e) 

l+Uc 


7 
7. 


(7c < 7e < 7m) 
-l/(l + at) 


[max( 7 c ,7 m ) < 7e] 
(16) 


In the limit at —>■ 0, one recovers dAfe/d 7 e oc 
[1 + y(7e)]“^7(r^ for 7c < 7e < 7m or dAfe/d 7 e oc 
[1 + y ( 7 e)]~^ 7 ~^~^ for 7 m < 7 c < 7 e; both match the ex¬ 
pressions derived in Nakar et al. (2009) and in Wang et al. (2010) 
for a homogeneous (non-decaying) turbulence. 

It is also instructive to show that the above average electron 
distribution, when associated with the appropriate (i.e. Lorentz fac¬ 
tor dependent) magnetic field, leads to the same scalings for the 
synchrotron spectrum as Eqs. l |12^ and l |13| >. Consider for instance 
the regime 7 > max ( 7 c, 7 m): the average synchrotron power 
of the blast scales as oc ( 7 s( 7 e) 7 e dA^e/dln 7 e; using 

the scaling for (7s ( 7 e) given in Eq. one recovers vF^^syn oc 


[1 -f y ( 7 e)] ^ 7e ^ indicated in Eq. 1 12 '. One can proceed simi 


larly for the other two regimes, noting in particular that for 7 m < 
7 e < 7 c (slow cooling), Ub = Ub- and it no longer depends on 

7c- 

Einally, using the above electron distribution, one can compute 
the inverse Compton component using standard formulae: 

-F^,ic(^'ic) Oi TIC j dq{l-u)gKi^{q) 

^ Fi/,svn -I 71 : T 0(1 tr) , 


47^5(1 — u) 


(17) 


with M = (1 -b z)hvicl (rb7emeC^); the function (;KN(q) = 
[^2q In q -1- (1 -t- 2q) (1 — 9 ) + 0^(1 — q)/ [2(1 -t- (7)]], with G = 
u/{l — u) characterizes the energy dependence of the Klein- 
Nishina cross-section, see Blumenthal & Gould (1970) for details. 

The prefactor ric defines the optical depth to Compton scat¬ 
tering; an explicit calculation leads to 


TIC = ^ 12crTrbncfdyn (18) 

Note that this expression uses A^e = cfdyn dTbn, i.e. it only 
includes the electrons that have been swept up in the last dynamical 


timescale, because electrons injected at earlier times have been adi- 
abatically cooled to Lorentz factors < 7 c and therefore do not par¬ 
ticipate in shaping the inverse Compton spectrum. With the above 
value of Tic, one can verify that the energy density of the radi¬ 
ation associated to the inverse Compton component is correctly 
normalized to yc times that contained in the synchrotron com¬ 
ponent, as can be verified by calculating J di'icTbic.ic in terms 
of f diZsynTj/Byn.syn (neglecting the influence of Klein-Nishina ef¬ 
fects). 


2.4 Slow cooling 

2.4.1 General procedure 


As discussed in detail in Nakar et al. (2009) and Wang et al. (2010), 
the cooling Lorentz factors and Compton parameters 7 c and yc can 
be obtained in the slow cooling regime 7 m < 7 c from the system 


yc (1 + yc) 


7c 


ee f ^ ^syn(< i>c) 

es- V 7 m J (7ayn(< Oa) 

7c,syn 

i-fyc 


(19) 


with 7 c,syn = (3/4)mec/(crT(7B-fdyn) the cooling Lorentz fac¬ 
tor in the relaxed turbulence, in the absence of inverse Compton 
losses. 

Write Cc = (7syn(< £'c)/(7syn(t'c) the term entering the first 
equation. If f-c < t^c, Cc < 1 because the peak of the synchrotron 
flux lies at or above (see below); in this limit, Klein-Nishina sup¬ 
pression inhibits the cooling of 7 c electrons. Assuming that such 
electrons cool by interacting with the Omin < o < Oc part of 
the synchrotron spectrum, which is generically the case, one can 
write Cc = oc 77 ^*'^”^^^^ and the system can be 

solved easily. 

In the opposite limit, i>c > Vc\ one may have Cc ~ 1 if the 
peak of the synchrotron flux is located at i/c, or Cc > 1 if the peak 
lies at higher frequencies. In this latter case, Cc — {vc/, 
where /3 is the synchrotron spectral index defined by Tj/.syn oc 
i/~^ in the spectral range above calculated thereafter. The ana¬ 
lytical calculation proposed here uses Eq. \21\ below to derive this 
P, but the result depends little on this choice, because in all cases 
considered, 1 — /3 is small, meaning that the peak flux is not very 
different from the flux vFi, at v^. 

The next step is to determine the critical Lorentz factors 7 c = 
7 (vc) and 7 m = 7 (i''m) with the general definition (Nakar et al. 
2009): 


7(iv) 


TbmeC^ 
(1 + z)hv 


( 20 ) 


which corresponds to the Lorentz factor for which electrons inter¬ 
act with photons of frequency u at the onset of the Klein-Nishina 
regime. Note that and are to be calculated in the relaxed 
turbulence of strength 5B- in this slow cooling regime. Using 
Eqs. ([^, ( | 10 ^ and (jTTJ, one may then calculate the Compton pa¬ 
rameters y (%) and y ( 7 m). Note also that the slow cooling limit 
7 m < 7 c implies % < 7 m. 

Another critical Lorentz factor is 70 , for which Y ( 70 ) = 1. If 
y ( 7 m) > 1 , then 7 m < 70 , and 70 can be obtained by solving 


y(7o)[i + y(7o)]~°*^^^+°*^ 

y(7m)[i-by(%0]““‘^‘'+“‘^ 


^( 70 ) 


-, 4/3 


7o 


at/G+at) 


which derives from Eq. {n}. This latter equation assumes that 
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i'( 7 o) lies in the spectral range between t/a (synchrotron self¬ 
absorption frequency) and Vm, in general a very good approxi¬ 
mation; it can be generalized to other cases without difficulty. If 
i^( 7 m) < 1 , then 70 < 7 in; one can repeat the above exercise 
to derive 70 , replacing the 4/3 exponents with (3 — p)/2, which 
characterizes the spectral dependence of between Vm and 

Vc- 

For improved accuracy, one may also determine the next order 
critical Lorentz factor 7 ^. = j [(1 -F z)hvsyn (%)]. Fol¬ 
lowing a procedure similar to the above, one can derive and 

The following makes use of the short-hand notation: v = 
Vsy-n ( 7 ) for various critical Lorentz factors; similarly, t'ayn( 7 o) is 
written j/q for commodity. 

One can derive the power-law 1 — jS' index of zzf/,syn as fol¬ 
lows. Below t/c, the scaling remains unchanged compared to the 
standard case because the cooling history does not influence the 
synchrotron spectram, and the results will not be repeated here. 
Above the spectral slope is determined by the scalings of Y ( 7 ), 
v{'y) and uFu^ayn oc (1 -F Y)~^ 7 ^”’’ [Efl- l|l2|]. In particular, if 
£'( 7 ) lies in a spectral domain in which i^F/,ayn oc , below 
the peak of the synchrotron energy flux, then electrons of Lorentz 
factor 7 cool by inverse Compton interactions with that portion of 
the synchrotron spectmm if Y ( 7 ) 3> 1, in which case Eqs. 

Go) and GD lead to 

Y oc 

2-8atl‘2 

V OC ^ 

fcool oc 7 “'^ ( 22 ) 

Assuming izFI,,ayn oc in the range of interest around iz( 7 ), 

and using the above scalings, one obtains directly 


2 - p - Qf -F (1 - /3)(1 -F at) 
2-at/2-F(l-^)Qt/2 


(23) 


The limit /3 —>■ 1 recovers the case in which v lies above the peak 
of the synchrotron flux, discussed previously. 


2.4.2 Power-law segments 


As mentioned above, the synchrotron spectrum at iz < i/^ re¬ 
mains unaffected and it is not discussed here. At the upper end 
of the spectral range, i.e. izo < v, Y ( 7 ) < 1 because of the 
Klein-Nishina suppression of electron cooling, which implies that 
t^Fl/.ayn OC with 


1-/3' 


2-p 

2 -at/[ 2 (l-Fat)] 


[vo < v\ 


(24) 


as can be derived from Eq. ( |12^ with 1 -F ^ — 1. This scaling also 
matches that derived by a full computation of the electron cooling 
history with negligible inverse Compton losses in Lemoine (2013). 
For at = 0, one recovers of course the standard fast cooling index 

P' = pI‘2. 

In the range izm < v < jzq, if it exists, the synchrotron 
spectrum is shaped by electrons with Lorentz factor 7 such that 
7 m < 7 < 70 , which thus cool by interacting with photons in the 
range v < Um\ one can therefore use Eq. ( |23f with 1 — /3 = 4/3, 
so that 


1-/3' 


5 - 3p/2 -F at/2 
3 -F at/4 


(zZm < V < Vo) 


(25) 


The limit at —> 0 gives P' = —2/3 -F p/2, which fits the results 
of Nakar et al. (2009) in this frequency range. 


In the range max (iZc, Pc) < v < min (Pm, tzo), the Lorentz 
factor of electrons shaping that part of the spectrum satisfies 
max( 7 c,%) < 7 < min ( 7 m, 70 ), hence one can use Eq. lj2^ 
with I — P = (3 — p)/ 2 , leading to 

l_ B' ~ 7- 3p-F (1 -p)at 

~ 4-F(1—p)at/2 

[max(zzc,Pc) < V < min (Pm, vo)] (26) 


Note that Pc < v^ is equivalent to Pc < v^. Here as well, one 
recovers the index j3' = —3/4 -F 3p/4 derived in Nakar et al. 
(2009) in the limit at —>■ 0. 

If Pc < Vc, meaning % < 7 c, the above completes the 
description of the spectrum. If, however, zzc < Pc, one needs to 
describe the intermediate range Vc < v < Pc. This range can 
actually be decomposed into two sub-ranges, as follows. 

For Pc < iz < Pc, the corresponding Lorentz factor verifies 
7 c < 7 < 7 c, hence iZc < P < Pc. Consequently, the particle 
cools on the spectral range of vF^^syn that it contributes to through 
synchrotron emission, so that one can use Eq. i23 1 with /3 = /?', 
giving 


l-;3' ~ 


— 2 -F 3at -F [4 -F (4 — 8p)at -F a?] 


211/2 


2at 


( 


Vc < V < Vf 


) 


(27) 


and 1 — P' 0.12 for at = —0.5, p = 2.3, i.e. a slowly rising 
uF 

^ i/tSyn- 

Finally, in the remaining range Uc < v < Uc, the particle 
cools on photons of frequency in the range Pc < P < Pc, for 
which 1 — p is given by Eq. l |26| l above. This leads to 

1 o/ c. 15 - 7p -F at(l - p) (10 - p -F at) /2 

P - 8 + at(l-p)(10 + at)/4 ^ ^ 

One finds 1 — ~ 0.09 for at = —0.5 andp = 2.3. 


2.4.3 Inverse Compton component 

In principle, one can derive analytical approximations to the inverse 
Compton component using the above broken power-law approxi¬ 
mations, folded over the particle distribution, as in Eq. GD above. 
However, this appears rather intricate, given the number of poten¬ 
tial power-law segments, in regards of the quality of the approxima¬ 
tion that one can obtain; indeed, as discussed in detail in Sari & Esin 
( 2001 ), folding over the distribution function generally introduces 
logarithmic departures, which smooth out the power-law segments 
and breaks. One can nevertheless describe the general features of 
this inverse Compton component in the slow cooling regime as fol¬ 
lows. 

Below tzic,m = 27miZm, the slope is that of vFi,^ayn below 
Vjn, i.e. 4/3. Eor vic.m < v < vic,c with vic,c = 2'ycVc, the 
slope of vFv.ic is (3 — p)/2, reflecting that of vF^^ayn in the corre¬ 
sponding range. Above izjc.c, the inverse Compton component re¬ 
flects, up to the afore-mentioned logarithmic corrections, the slope 
of vFv.ayc above Vc, which is generally close to flat or gently ris¬ 
ing, see above. 

Consequently, if the synchrotron spectrum has an extended 
range above Vc where it is close to flat (i.e. /3 ~ 1), then one 
might have a close to flat inverse Compton component, at least up 
to the cut-off frequency defined as 

i2lC,KN = TcT'c 


(29) 
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corresponding the boosting of photons at the onset of the Klein- 
Nishina regime by 7 c electrons. One can define another cut-off fre¬ 
quency, as follows: 

= 7ct^c (30) 

which corresponds to the boosting of photons by electrons of 
Lorentz factor 7 c, at the onset of the Klein-Nishina regime. The 
ratio t'lc.c/^'ic.KN can be written as 7c/7c, hence the ordering of 
one with respect to the other depends on whether 7 < 7 c or not. 
In any case, the actual cut-off occurs at j^ic.kn, while the pres¬ 
ence of r'ic.e may lead to a feature (e.g. softening) in the inverse 
Compton spectrum. This can be understood by noting that in the 
present case, the synchrotron (energy) flux generically peaks above 
Vc, while the particle distribution function falls steeply beyond 7 c, 
hence the peak of the inverse Compton component is determined 
by the boosting of photons by electrons of Lorentz factor 7 c. 


2.4.4 Comparison to numerical calculations 


The above analytical broken power-law model of the synchrotron 
spectrum is compared to a full numerical calculation (with the algo¬ 
rithm described in Sec. |2.5| thereafter) in Fig.[^ in two different rep¬ 
resentative cases: upper panel, observer time fobs = 10'* s, blast 
energy E = 10®® ergs, external density n = 0.01 cm“®; lower 
panel, fobs = 3 x 10* s, E = 10®* ergs, n = 10®® r~^ cm“® 
(wind profile with shock radius r expressed in cm); for both, 
Ce = 0.1, p = 2.3, tB+ = 0.01 and at = —0.4, assuming 
cb = CB+ [f/(100a;“;*)] . For a decelerating adiabatic Bland- 

ford & McKee (1976) solution, the value of the blast Lorentz fac¬ 
tor at these observer times are Fb — 33 for the upper panel and 
Fb — 29 in the lower panel. For the above decay law of the mag¬ 
netic field, one finds for the first scenario cb- = 3.2 x 10“® 
(fdyn = 1.3 X 10® s, cUpj* = 7.6 X 10“® s), and in the sec¬ 
ond scenario e_B_ = 2.1 X 10“® (fdyn = 1.7 X 10® S, CJpj* = 
3.5 X 10"® s). 

The critical frequencies are indicated with arrows. The thick 
solid line corresponds to the numerical calculation (synchrotron in 
blue, inverse Compton component in orange) while the dashed line 
shows the analytical estimates, which clearly provides a faithful 
match in both cases. 

In the first scenario (upper panel), % ~ 60 while 7 c ~ 10^: 
Klein-Nishina effects are therefore particularly strong; there are ac¬ 
tually so strong that j/q < r'c, which means that the Klein-Nishina 
suppression of electron cooling reduces the Compton parameter 
to below unity at 7 c. Consequently, the dependence of K on 7 
does not affect the spectrum above v^. Noting that in this region 
vFi, oc 7 ®"^ [Eq. 1 12 ] and u oc 72 -“t/ 2 (i+“t) |-gq ([7o|], one 
finds 1 /FI, oc 1 /*“^ with 1 —/3 = (2 — p)/[2 — at/2(l-I-at)], 
which matches Eq. \2A) . As mentioned earlier, the spectrum re¬ 
mains unaffected with respect to the standard synchrotron spectrum 
below Vc because the electrons shaping that part of the spectrum do 
not cool on a dynamical timescale. 

In this first scenario, the method proposed in Eqs. fig over¬ 
estimates 7 c by a factor 2.5, hence by a factor 6 and t'ic.KN = 
Fc7c by a factor 2.5 as well. Taking into account this overestimate, 
the numerical calculation indicates that the suppression of the in¬ 
verse Compton flux becomes noticeable at a factor ~ 5 below the 
theoretical value of zzic.kn, as calculated with the correct 7 c. As 
anticipated, the characteristic frequency I'ic.e leads to a soft soft¬ 
ening in the inverse Compton component, but not to a cut-off. 

In the second scenario (lower panel), 7 c < % hence Klein- 


Nishina suppression of the inverse Compton cooling becomes ef¬ 
fective at Pc only. The (analytical) synchrotron spectrum is thus 
close to flat in the region iZc i' ^ Pc, as indicated by Eqs. |28]l, 
above, then rising with 1 — /? — 0.14 corresponding to 
Eq. l |26) above for the range Uc < v < uq, because ug < Om. 
Above r/Q, Eq. \2A\ applies and gives the same high energy spectral 
slope as for the previous scenario. 

In this second case, the analytical calculations underestimate 
7 c by a factor 1.8. There is nevertheless a broad satisfactory agree¬ 
ment between the analytical synchrotron spectrum and the numeri¬ 
cal calculation. 


2.4.5 Comparison to non-decaying scenarios 


Figure 1^ provides a numerical comparison of the spectra shown 
in the lower panel of Fig. with two calculations for the same 
parameters but a homogeneous (non-decaying) turbulence: one in 
which SB ~ CB+ = 0.01, another one in which cb = es- = 
2.1 X 10~®, which corresponds to the value of es(fdyn), i.e. close 
to the contact discontinuity, in the above decaying micro-turbulence 
model. 

As expected the SSC spectrum with decaying microturbu¬ 
lence merges with that corresponding to uniform cb- at frequen¬ 
cies below v^, since electrons of Lorentz factor 7 < 7c then cool 
in magnetized turbulences of equal strength in both models. The 
synchrotron spectrum for decaying micro-turbulence also merges 
with the synchrotron spectrum for uniform tB+ at the highest fre¬ 
quencies, since the cooling time for those emitting electrons be¬ 
comes shorter than A, hence the particles effectively cool in a 
magnetic field characterized by tB+. However, the integrated pow¬ 
ers for these two models differ, because the cooling efficiencies 
~ (7c/7m)^“^ differ. 

In this regard, the slow cooling synchrotron spectrum for de¬ 
caying micro-turbulence is a hybrid of the spectra for uniform high 
and low es, transiting from eb- at low frequencies to £ 3 + at high 
frequencies. This justifies the use of a two zone model, one with 
low Efl- and one for high £b-i-, to compute an approximated spec¬ 
trum in wavebands at respectively low and high frequencies. 

In Fig. the synchrotron spectra have been arbitrarily con¬ 
tinued at very high frequencies, albeit with a thin line, but they 
should of course cut-off at some maximal energy where the acceler¬ 
ation timescale becomes of the same order as the cooling timescale. 
Since this depends on acceleration physics, see the discussion in 
Lemoine (2013), the lines have been turned from thick to thin at 
an ad-hoc location corresponding to a synchrotron photon energy 
of 1 GeV. This estimate is discussed in Plotnikov et al. (2013), 
Lemoine (2013), Wang et al. (2013), Sironi et al. (2013); it depends 
on the afterglow parameters and, in particular, on observer time. 

For reference, one notes the critical frequencies: 


izic.c 


I'IC.KN 


27®iZc 

7.3 X 10®®Hzii;54£B!.(l5^:,n'7iob=,4.5^;'3"c:2" 
1 P 2 

——Fb 7 cmeC 

14-2 


— 4 X iU nz 41/54 CB-,-5^*,11.7^ohB,4.5^-t- ^c,2 

(31) 

with the notations z+ = (1 -F z)/2, 17,2 = (1 -F 17)/100 and 
/4 *,ii.7 = A*/(5 X 10** g/cm®). Forreference, in the scenario of 
Figs.||and[^ A* ~ 0.3 and Kc 50fores- = 2.1 x 10"®. Note 
the strong dependence of vic,c and vkn on the external density. 


3/2 ,1/2 


,-3/2v-l 
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Figure 1. Comparison of the analytical calculation (dashed red line) of the synchrotron spectrum in the slow cooling regime to a numerical calculation (solid 

blue line), for two representative cases: upper panel, observer time fobs = 10^ s, blast energy = 10^^ ergs, external density n = 0.01 cm~^; lower panel, 
= 3x10'^ s,F? = 10^^ergs,n = 10^^ cm“^; in both cases, eg = 0.1,p = 2.3and€s = ^b+ . The analytical estimates 

of the characteristic frequencies are indicated with arrows. The solid orange line represents the numerical calculation of the inverse Compton component. 



Figure 2. Synchrotron and inverse Compton spectra at fobs = 3 x 10“^ s, representative of the slow cooling regime, for a blast with energy E = 10^^ ergs 
impinging on a progenitor wind with density n = 10^^ r~^ cm~^ (r in cm), assuming eg = 0.1 and p = 2.3, for three microphysical models, as 

indicated: homogeneous (non-decaying) = ^b+ = 0-01 (dashed red line); decaying = ^b+ (solid blue); homogeneous 

eB = es- = 2.1 X 10“^ (dash-dotted orange), the value of es- being representative of es close to the contact discontinuity in the decaying es model. 
The synchrotron model predictions have been thinned beyond an ad-hoc maximal synchrotron photon energy of 1 GeV (see text). Characteristic frequencies 

are: for at =—0.4 and for eg = — 5.6 x 10^^ Hz, i/c — 1.2 x 10^® Hz and t/iC,KN — 2 x 10^® Hz; for = es-i-, ^'m — 1.2x10^^ Hz, 

i/c ~ 4.2 X lO^^Hzandz/ic.KN ~ 2 x lO^^Hz. 
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As discussed above, the peak of the inverse Compton component 
is expected to occur at !^ic,kn, although the numerical calculation 
suggests that the turn-over becomes manifest a factor ~ 5 below 
the above theoretical value. 

The VERITAS collaboration has recently heen able to observe 
the exceptional GRB130427A and to put stringent upper limits on 
the emission at > 100 GeV (Aliu et al. 2014). The absence of de¬ 
tection of this energy range suggests that, for this burst at least, the 
inverse Compton component has cut-off below ~ 100 GeV, while 
the Fermi detection of multi-GeV photons up to a day or so sug¬ 
gests that this cut-off lied above 1 — 10 GeV. Such a cut-off energy 
fits well with the above estimates for z/ic,kn for a low average cb. 

2.5 Fast cooling 

The fast cooling regime involves a susbtantial variety of syn¬ 
chrotron spectra, with multiple breaks and indices, as discussed in 
detail in Nakar et al. (2009) and Wang et al. (2010) for the case 
of uniform cb • One key difference with the slow-cooling regime is 
that particles with Lorentz factors 7 < max ( 7 c, 7 m) may have a 
non-trivial cooling history while in the slow-cooling regime, such 
particles do not cool. As a consequence, it is difficult to even derive 
the cooling Lorentz factor and the Compton parameter Yc when 
Klein-Nishina effects become significant. 

In this fast-cooling regime, it is actually more efficient to com¬ 
pute the spectrum numerically, using the following simple and ef¬ 
ficient algorithm. One starts with a template synchrotron spectrum, 
for instance that corresponding to a homogeneous magnetized tur¬ 
bulence. One can then derive a first approximation to 7 c and Yc, 
either using standard formulae (e.g Panaitescu & Kumar 2000) - 
which ignore KN effects - or through an explicit determination of 
7 c as the Lorentz factor for which cooling takes place on a dy¬ 
namical timescale, using a radiation energy density inferred from 
the template spectrum. With 7 c and Yc, one can solve Eqs. 

and (n} to compute the frequencies and Compton parame¬ 
ter as a function of the initial Lorentz factor of an electron; one 
can then use Eqs. © and to compute an improved version 
of the synchrotron spectrum, properly taking into account the cool¬ 
ing history of the electrons in the decaying turbulence as well as 
all relevant Klein-Nishina effects. This latter spectrum remains an 
approximation, because it relies on a guessed value for jc and Yc. 
Nevertheless, iterating the above process, using each time as a tem¬ 
plate the previously computed synchrotron spectrum, one obtains 
after ~ 10 iterations a self-consistent synchrotron spectrum, with 
7 c and Yc determined to high accuracy. Finally, one can derive the 
inverse Compton spectrum using Eq. l |17^ . 

The above algorithm provides a self-consistent estimate of the 
synchrotron and inverse Compton spectra with a normalization ac¬ 
curacy of order unity. This accuracy can be checked by calculating 
a posteriori the integrated synchrotron and inverse Compton pow¬ 
ers and comparing to the total electron power injected through the 
shock; in the fast cooling regime, these should match. The error 
is of order 10 — 40 % for the SSC spectrum of a decaying micro¬ 
turbulence as shown in Fig.[^ depending on observer time; it is less 
than 10 — 20 % for at = OandcB = cb-i-, but it becomes a factor 
<2 for at = OandcB = cb--M ost of the error results from the 
broken power-law normalization of the flux in Eq. \\2\ and from 
the treatment of the Klein-Nishina cross-section as a step function 
in the calculation of the synchrotron cooling history. The resulting 
uncertainty remains nevertheless satisfactory given the uncertainty 
associated for instance to the definition of t^yu (hence 7 c) in the 
absence of a realistic description of the blast energy profile. 


A detailed example of the SSC spectrum of the blast, for the 
same parameters as in Fig.|^ is presented in Fig. at an observer 
time fobs = 30 s; assuming a Blandford & McKee (1976) decel¬ 
erating solution in a wind profile, the Lorentz factor of the blast at 
that time is Tb — 160. The value of cb- in this case is 4.2 x 10“®. 
As expected, the spectrum corresponding to a decaying turbulence 
merges with the spectrum for uniform es = £ 3 + above a fre¬ 
quency u ~ 10^^ Hz, since the electrons that emit in that range 
cool fast, in a region where cb — eB+. In this fast cooling regime, 
the total integrated energy densities of the three SSC spectra cor¬ 
respond to the injected electron energy density. The spectrum for 
decaying micro-turbulence does not merge with that for uniform 
cb = cb- at low frequencies, since the cooling Lorentz factors 
differ for both. In this fast cooling regime, one cannot therefore de¬ 
scribe accurately the synchrotron spectrum at low frequencies with 
a spectrum computed for uniform low cb : an explicit calculation 
becomes necessary. 

Finally, note that the inverse Compton component in the fast 
cooling regime is expected to peak at 'y^Um, if the synchrotron flux 
peaks at Vm and if one omits KN effects. The Klein-Nishina sup¬ 
pression implies a turn-over of the IC component at most at 7 mf'm, 
for reasons analog to those discussed in the slow cooling regime, 
see also Nakar et al. (2009). If Vui < Tm, the slow rise of the 
synchrotron flux above in the case of low cb (due to the KN 
suppression of electron cooling) or decaying microturbulence im¬ 
plies a comparable behavior of the IC component between 7 mi^m 
and 7 mTm- This feature is not clearly seen in Fig.j^due to the (rel¬ 
ative) proximity of these two frequencies, 7^i>m 6 X 10^® Hz 

and7mi/in ~ lO^'^Hz. 


3 SPECTRA AND LIGHT CURVES 
3.1 Spectral and temporal behaviors 

The spectra shown in Figs. and [^illustrate how a complete and 
simultaneous spectral coverage would allow to tomograph the evo¬ 
lution of the micro-turbulence behind the relativistic shock. The 
effects are most noticeable in the X-ray and MeV regions, as one 
would expect: in this region of the spectrum, the emitting electrons 
feel the decaying turbulence, while at the highest frequencies, they 
cool in regions with e ~ £b+, and at frequencies v < Vc they 
cool in a magnetic field characterized by eb-, the value of which 
evolves slowly in time. 

Afterglow models often rely on the spectral and temporal 
slopes in various domains and their so-called closure relations to 
make comparison to observations. Figure [^therefore presents the 
spectral slopes /? defined by oc (where Fi, sums the syn¬ 
chrotron and inverse Compton fluxes), for the optical. X-ray, high 
energy (0.1 — 10 GeV) and very high energy (> 10 GeV) ranges 
(note that no attenuation in extra-galactic background radiation has 
been assumed for the latter range). This figure compares the spec¬ 
tral slopes for the three previous representative models: (at = 0 , 
£b = 0.01), at = —0.4 and (at = 0, eb = 10“®) with other¬ 
wise same parameters as in Figs. min except k, which takes values 
0 (constant density profile) or 2 (stellar wind). 

In the optical, {3 has been calculated at a reference frequency 
of 4.7 X lO^'^ Hz (R band); in the X-ray, /3 is calculated as the 
average slope over the interval of energies 0.3 — lOkeV; at high 
energy, it is calculated as the average of the energy interval 0.1 — 
10 GeV and at very high energy, over > 10 GeV. 

According to Fig. [^ the most robust signature of a decay¬ 
ing micro-turbulence appears to be a slightly harder slope in the 
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Figure 3. Same as Fig. I^at observer time tobg = 30 s, representative of the fast cooling regime. Characteristic frequencies are: for = —0.4,1^0 ~ 10^^ Hz, 
t'm ~ 3 X 10^® Hz and 7^!Zm ~ 6 X 10^® Hz; for es = cs—, r'c — 2 X 10^® Hz, Vui ~ 3 X lO'^® Hz and 7^f'm ~ 6 X 10^® Hz; for eb = es+, 
t/c ~ 4 X 10“ Hz, i/m ~ 4 X 10“ Hz and 7^£;in ~ 6 X lO^® Hz. 


X-ray, P ~ 0.9 [panels (b) and (e)] vs /3 ~ 1.15 for uniform 
e_B in the first hours [panels (a), (c), (d) and (f)]. The latter value 
corresponds to the fast cooling regime /? = p/2, therefore it de¬ 
pends on p; however, one does not expect it to go below 1, because 
p > 2 is a generic prediction of relativistic shock acceleration, e.g. 
Bednarz & Ostrowski (1998), Kirk et al. (2000), Achterberg et al. 
(2001), Lemoine & Pelletier (2003), Sironi et al. (2013). Current 
data do not allow to distinguish between these limits; in particular, 
the Swift data lead to /3 ~ 1 ±0.1 (Evans et al. 2009) in afterglows 
with standard power-law decay. Interestingly, even in the case of a 
homogeneous turbulence, the X-ray slope hardens at late times be¬ 
cause of the emergence in the X-ray range of the inverse Compton 
component; for instance, in panels (a) and (c), one can see /3 transit 
to values of order 0.7, corresponding to the low energy extension 
of this inverse Compton component with slope /3 = (p—l)/2. 

In the optical range, the slope is comparable to that in the X- 
ray range for the decaying micro-turbulence scenario at observer 
times ~ 10® s, but significantly harder at earlier times when the 
optical falls in the range i/c — iZm: Fig. [^indicates values/? ~ 0.3, 
harder than expected (1/2) in the standard fast cooling regime in 
this range of frequencies. At an observer time fobs = 30 s, for 
fc = 2 corresponding to panel (e) as well as to Fig.[^ — 1.4 x 
10®® Hz, which lies below the peak of the synchrotron component 
(see Fig.[^. This implies that Klein-Nishina effects are significant, 
and that 7 ni electrons cool by interacting with the segment in the 
range 10®^ — 10®® Hz, whose index ~ 0.8. In this case, one can 
compute the expected index of the segment below zzm, using a 


variant of Eq. ([2^: one notes that izF„_syn oc 7 /(l-|-F) [Eq. 1 13 i] 
which gives 


1-/3' = 


1 — Qt -f (1 — /9)(1 -b at) 

2-atl2 + {l- P)atl2 


~ 0.3 


(32) 


the last equality applying for at = —0.4 and /? = 0.8. This 
explains the values of at found in the optical range. Such values 
do depart from the standard synchrotron spectra, although Klein- 
Nishina effects may also cause values /3 ~ 0.3 below Vuiin in the 
case of homogeneous turbulence scenarios, see Nakar et al. (2009), 
their figures 2 and 3 for example. Therefore, it is not clear at present 


whether one can consider such values of /3 as a clear signature of a 
decaying micro-turbulence. 

All in all, an accurate measurement of /3 in the X-ray range, or 
better, in the MeV range if the MeV afterglow could be detected, 
would provide the best probe of at, see also Fig.|^and[^for illus¬ 
trations of these effects. 

Another quantity of interest for a general description of the 
afterglow is the temporal slope, defined by T)/ oc Of course, 
these temporal slopes directly depends on the time evolution of the 
various parameters, through the assumed evolutionary law for Tb 
and the evolution of r and n, in contrary to the spectral slopes /? at 
any given time. Here Tb is assumed to decrease as in the Blandford 
& McKee (1976) adiabatic solution, i.e. Tb oc 

The values of a for the same models and intervals as Fig.|^are 
reported in Fig.|^ The largest differences between constant (high) 
and decaying es result from the different transit times between the 
slow and fast cooling regimes, but these times depend in turn on 
other parameters that are a priori unknown. The light curves other¬ 
wise present similar features, without a clear trend distinguishing 
one from the other. 

Figures 12 and thus indicate that a and P are by themselves 
weakly sensitive probes of the dynamics of the magnetized turbu¬ 
lence in the blast and that it is not possible at present to distinguish 
a decaying micro-turbulence from a uniform low- or high es on the 
basis of these data. It appears much more effective to try to probe 
es through a multiwavelength fit of the afterglow light curves, us¬ 
ing not only the temporal and spectral slopes, but also the ratio of 
fluxes between different spectral windows, as done in Lemoine et 
al. (2013), Liu et al. (2013) for Fermi-LAT bursts. 


3.2 Emission at very high energy 

Finally, an interesting consequence of a decaying micro-turbulence 
is the generic prediction of substantial emission at the highest en¬ 
ergies, due to the large value of the Compton parameter. Naive 
estimates, Y ~ indicate values of several hundreds 

for y, although they neglect Klein-Nishina effects which depend 
on the electron energy, therefore on observed frequency; further- 
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Figure 4. Flux density index /3, with Fi, oc , as a function of time, 
in various wavebands, for different external density profiles: k = 2 (wind, 
n = cm“^) and k = 0 (constant density, n = 1 cm“®), assum¬ 

ing uniform eg = 0.01, es = (100ajpit)“‘ or uniform eg = 10“®. 
Blast energy, jet Lorentz factor and Ce are as in Fig.l^ Wavebands are as in¬ 
dicated: for reference, optical (circles) u = 4.7 X 1(F'^ Hz, X-ray (squares) 
integrated from 0.3 to 10 keV (i.e. 0.72 — 24 X 10^^ Hz), 0.1-10 GeV (up¬ 
ward triangles, in frequency 0.24 — 24 X 10^® Hz) and > 10 GeV (dia¬ 
monds, in frequency > 24 X 10^^ Hz). 
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Figure 5. Observer time decay index of the flux density, i.e. Fi, oc 

for various models, as indicated in Fig. Blast energy and Ce are as in 
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Figure 6. Ratio of the total energy flux above 30 GeV to a theoretical flux 
obtained by matching the calculated total flux at 0.1 GeV and extrapolating 
this flux to higher energies with a power-law ^ , as in Inoue 

et al. (2013). The symbols correspond to the six models studied in Figs.j^ 
and[^ as follows: red circles €b = and at = 0 ; blue squares 

at = —0.4; orange triangles eg = eg_ and at = 0. Filled symbols 
con'espond to fc = 2 and open symbols to k = 0 ; other parameters are as 
in Fig.[^ 


more, the non-trivial spectral shape above t/c in the case of decay¬ 
ing micro-turbulence modifies the ratio of the inverse Compton to 
the synchrotron component. As discussed in Wang et al. (2013), 
emission above 10 GeV is most likely of inverse Compton origin, 
because the maximal synchrotron photon energy is more likely of 
the order of 1 GeV or so at 100 — 1000 sec observer time. This 
emission is a prime target for future gamma-ray telescopes such as 
HAWK (Mostafa 2013) or CTA (Inoue et al. 2013). For the partic¬ 
ular case of CTA, Inoue et al. (2013) have investigated the detec¬ 
tion rates of gamma-ray bursts above 30 GeV by assuming that the 
spectrum continues beyond IGeV with a spectral index /3 = 1.1 
(corresponding to the standard fast cooling regime /3 = p/2 with 
p = 2.2) and scaling the flux at 1 GeV to that measured by the 
Fermi-LAT instruments. Their simulations lead to about one de¬ 
tection per year. This rate is rather low, therefore any improvement 
would be quite valuable, given the potential impact of a high energy 
detection. 

One can use the calculations of Sec.|^to study how a decay¬ 
ing micro-turbulence or low es_ affects these predictions. In or¬ 
der to do so, one calculates the ratio of the total (synchrotron -l- 
inverse Compton) energy flux above 30 GeV, F{> 30 GeV) = 
/soGeV ® theoretical reference flux Tth(> 30 GeV). As 

in Inoue et al. (2013), Fth(> 30 GeV) is obtained by extrapolating 
the total flux measured at a reference energy, here 0.1 GeV, with an 
index/3 ~ 1.1, i.e. 

Fth(>30GeV) = F^O.lGeV) 

xT d4^^V"(33) 

J 30 G.V/H \0.1GeYj 

with Fv = Fv^syn + Fv,ic the total flux. The reference energy 
chosen here is smaller than that in Inoue et al. (2013), because the 
inverse Compton flux is already prominent at 1 GeV in the scenar¬ 
ios studied, as shown in Fig.|^and[^for example. However, given 
the value of the index /3, the theoretical flux vFv is roughly flat 
above 1 GeV, hence this should not affect the statistics of detection. 

The results are shown in Fig.[^ which carries out this evalua¬ 
tion for the six models shown in Figs.|^an d|5] This figure indicates 


that a decaying micro-turbulence with at = —0.4 (or a low aver¬ 
age Efl) increases by a factor of a few, up to an order of magnitude, 
depending on at, k and toha, the prospects of observing the after¬ 
glows at energies > 30 GeV, relatively to statistics computed for a 
model with es = 0.01, as in Inoue et al. (2013). This certainly 
brings the number of potential detections by instruments such as 
CTA in a more confortable range. 

Furthermore, it is important to note that in most models stud¬ 
ied here, the inverse Compton component contributes to a signifi¬ 
cant fraction of the total flux at 0.1 GeV, therefore the above ratio 
actually is an underestimate of the ratio of the inverse Compton flux 
at high energies to the synchrotron flux at GeV energies. 


4 CONCLUSIONS 

This paper has discussed the spectral shapes of the synchrotron- 
self-Compton spectrum of a relativistic blast wave, including all 
relevant Klein-Nishina effects, and their evolution in time. A par¬ 
ticular emphasis has been put on the impact of a decaying micro¬ 
turbulence behind the shock front, which is motivated by theoretical 
analysis (Chang et al. 2008, Lemoine 2015) and observational in¬ 
ference (Lemoine et al. 2013, Liu et al. 2013). However, the results 
are fully applicable to the case of a uniformly magnetized blast, 
possibly with a low value of the average es. 

A decaying micro-turbulence and/or a low average value of 
es both lead to a large Compton parameter, with interesting 
physical and observational consequences. From a more theoreti¬ 
cal point of view, the Klein-Nishina suppression of inverse Comp¬ 
ton cooling has a strong effect on the synchrotron spectral shape, 
as noted elsewhere for the case of a uniformly magnetized blast 
(Nakar et al. 2009, Wang et al. 2010). In the case of a decaying 
micro-turbulence, the modification is not trivial to compute, and 
the present paper has described a simple algorithm which allows 
to compute the full SSC spectrum with satisfactory accuracy, at a 
modest numerical cost. Among the interesting phenomenological 
consequences, one may point out the slight deviations in the spec¬ 
tral and temporal slopes induced by the decaying micro-turbulence, 
or by KN effects in uniformly magnetized blasts at low es. A multi¬ 
wavelength coverage of the afterglow would, in principle, allow one 
to tomograph the dynamics of this magnetized turbulence, through 
its influence on the light curves in various wavebands. However, 
as expressed in terms of the spectral /) and temporal a slopes, de¬ 
fined customarily through oc the deviations are rela¬ 

tively weak and not currently distinguishable through observations. 
A multiwavelength fit of the afterglow, which also relies on the flux 
ratios between various wavebands, seems to provide a more sensi¬ 
tive probe of the dynamics of the micro-turbulence. 

Finally, a large parameter also implies a large inverse 
Compton flux at multi-GeV energies, relatively to the lower energy 
synchrotron flux, with direct consequences for the detectability of 
gamma-ray bursts afterglows by future gamma-ray telescopes. A 
numerical estimate indicates that a low average es would imply a 
detection rate several times larger than currently anticipated on the 
basis of the extrapolation of the flux of gamma-ray bursts detected 
by Fermi-LAT. 
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